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Oj. We report the discovery of a tight substellar companion to the young solar analog 
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PZ Tel, a member of the p Pic moving group observed with high contrast adaptive 
optics imaging as part of the Gemini NICI Planet-Finding Campaign. The companion 
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was detected at a projected separation of 16.4±1.0 AU (0.33±0.01") in April 2009. 
Second-epoch observations in May 2010 demonstrate that the companion is physically 
associated and shows significant orbital motion. Monte Carlo modeling constrains the 
orbit of PZ Tel B to eccentricities >0.6. The near-IR colors of PZ Tel B indicate a 
spectral type of M7±2 and thus this object will be a new benchmark companion for 
studies of ultracool, low-gravity photospheres. Adopting an age of 12]']4 Myr for the 
system, we estimate a mass of 36ib6 Mj^p based on the Lyon/DUSTY evolutionary 
models. PZ Tel B is one of few young substellar companions directly imaged at orbital 
separations similar to those of giant planets in our own solar system. Additionally, 
the primary star PZ Tel A shows a 70 /im emission excess, evidence for a significant 
quantity of circumstellar dust that has not been disrupted by the orbital motion of the 
companion. 

Subject headings: brown dwarfs - planetary systems - stars: pre-main sequence 



1. Introduction 
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2003 ). largely at separations >50 AU. Only a handful of 



substellar or planetary mass companions have been found at separ ations <20 AU, including the 
recently confirmed planet imaged around /3 Pic (ILagrange et al.ll2010l ') and a few substel lar compan- 



ions around older (;^Gyr) stars (e.g., HR 7672 B [Liu et al.ll2002l |: SCR 1845-6357AB [Biller et al 
200fi| ). 



The Near-Infrared Coronagraphic Imager (NICI) at the 8.1-m Gemini-South Telescope (jChun et al 



20081 ) is a dedicated AO instrument tailored expressly for direct imaging of brown dwarf and exo- 
planet companions. NICI combines several techniques to attenuate starlight and suppress speckles 
for direct detection of faint companions to bri ght stars: (1) Lyot coronagraphy , (2) dual-channel 



imag i ng for Spectral D ifferential Imaging (SDI: iRacine et al.lll999l : iMarois et al.ll2005l : IClose et al 



20051 : iBiller et al.ll2007l) and ( 3) operation in a fixed Cassegrain rotator mode for Angular D ifferen- 
tial Imaging (ADI; Eiulbooi iMarois et aPboOfil : iLafreniere et al.ll2007bl : biller et ahlboosl). While 
these techniques have been used individually in large planet-finding surveys (jBiller et al.l 120071 : 
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Lafreniere et al.ll2007al : iNielsen k Closdboid ). the NICI Campaign is the first time ADI and SDI 
have been employed simultaneously in a large survey. 



Since December 2008, the NICI Planet-Finding Campaign (JLiu et al.ll2009l ) has been obtaining 
deep, high-contrast adaptive optics (AO) imaging of a carefully selected sample of ~300 young, 
nearby stars. We report the discovery of a close substellar companion to PZ Tel, a young solar 
analo g (KO spectral type) and a member of the ^12 Myr Pic rn oving group (jZuckerman et al 
200 ih with a Hipparcos distance of 51.5±2.6 pc (jvan LeeuwenI 120071 ). 



2. Observations and Data Reduction 

The NICI Campaign uses specialized observing strategies to exploit the unique capabilities 
of the NICI instrument. While SDI provides significant contrast gains even for non-methanated 
objects at separations down to 0.7", at ~0.3" PZ Tel B is significantly self-subtracted in our SDI 
data reductions. Thus, in this letter we focus solely on ADI-based techniques. For ADI, the image 
rotator is turned off, and the field is allowed to rotate with the parallactic angle on the sky. The 
image of a true companion will rotate with the parallactic angle, while speckles will stay fixed 
relative to the detector. A PSF for the primary star can be built and subtracted from each frame 
to remove the speckle pattern. 



2.1. First and Second Epoch ADI -|- SDI Observations 

We obtained ADI and SDI observations of PZ Tel on 11 April 2009 (UT) and 9 May 2010 (UT). 
These observations were conducted using a flat-topped Gaussian profile focal plane mask with 
HWHM=0.32" (referred to hereafter as the 0.32" mask), the CH4 4% Short filter {Xcentrai = 
1.578 fim) in the blue channel, and the CH4 4% Long filter {Xcentrai = 1.652 /xm) in the red 
channel. At each epoch, we acquired 45 1-minute frames in a combined ADI -|- SDI mode, over 
16.7° of field rotation in epoch 1 and 23.2° in epoch 2. 

Data for each camera were dark-subtracted, fiat-fielded, and distortion-corrected. Sky back- 
ground was negligible compared to halo brightness at the companion location and has not been 
removed. The focal plane mask is not completely opaque at its center, leaving an attenuated, un- 
saturated image of the primary star (henceforth "starspot") which we use for image registration. 
Images from both channels were registered to the first blue channel image according to starspot 
centroid positions. To increase our sensitivity to non-methanated companions, we performed a 
"broadband ADI-only" reduction by combining red and blue channel images instead of subtracting 
them. The median combination of all the images was subtracted from each individual image to 
remove the light from the primary star. Individual images were then rotated and stacked. 



2.2. 0.22" Mask JHKs Observations 

NICI coronagraphic observations using the narrower flat-topped Gaussian focal plane mask 
with HWHM=0.22" (hereafter 0.22" mask) were acquired in dual-channel mode at iJ-band and 
K^-band on 9 May 2010 UT and in single channel mode at J-band on 10 May 2010 UT The 



NICI filters are on the M auna Kea Observatories photometric system (ISimons &: Tokunagall2002l : 



Tokunaga Sz Vaccal l2005l ) . We acquired 10 60-s frames in each filter with the rotator on. Frames 
were dark-subtracted, flat-fielded, and stacked according to starspot centroid position. The DC level 
calculated in an annulus far from the object was subtracted to remove sky background. JHKs-hand 
images are presented in Fig. [TJ 

The 0.22" mask has an appproximately Gaussian transmission profile, an inner flat plateau 
of 0.12" radius, and an outer cutoff of 0.4" radius. With its separation of ~0.3", PZ Tel B is 
lightly attenuated by the mask. Thus we must determine the mask transmission at both the center 
and at the companion's position in order to measure the companion's flux ratio relative to the 
primary. To measure the transmission as a function of radius, we acquired JHKs-hand images 
of a 2MASS-selected 4.4" binary star. For each filter, the primary component was scanned across 
the mask along the same angular trajectory as PZ Tel B. Scans were begun with the primary 
under the starspot in order to accurately determine the mask center position. Images were taken 
at 0.02"steps out to 0.5". Integration times were chosen to keep both components unsaturated. 
We measured the flux ratio of the primary to the secondary as a function of the primary star 
displacement from the mask center using photometric apertures of 1 to 3 pixels in radius. For 
each photometry aperture, we fit a polynomial power-law to the mask transmission as a function of 
radius assuming a perfectly transparent mask beyond 0.43" radius. We used this fit to correct the 
measured flux ratio of PZ Tel B. At a separation of 0.33", we estimate an error in our calibration 
of the mask transmission of 6%, 3%, and 11% in J, H, and Ks respectively. These mask calibration 
errors are included in the final photometry errors. 



2.3. Maskless Narrow-band ADI Observations 

To acquire photometry and astrometry without attenuation due to the focal plane mask, we 
acquired 43 maskless narrow-band ADI images on 10 May 2010 UT. Each image had an integration 
time of 30.4 s, and the field rotated by 25.7° over the entire observation. The narrow-band CH4 
1% Short filter {Xcentral = 1.587 //m) was used in the red channel and the H2 2-1 narrow-band filter 
{^central = 2.1239 //m, width of 1.23%) was used in the blue channel. Data from each filter were 
reduced in separate ADI reductions (§ 2.1). 



To correct losses due to self-subtraction from ADI (e.g. iMarois et al.ll2006l ) and to determine 
the photometric and astrometric uncertainties, we inserted simulated objects into each individual 
frame by scaling a cutout of the stellar peak and shifting it into position. Objects were simulated 
at a separation of 0.36", at position angles of 0, 180, and 270° (counterclockwise from north), and 



with Amag = 4.6-5.8 in steps of 0.1 mag. Objects were retrieved using apertures of 1-3 pixels in 
radius to determine self-subtraction as a function of magnitude. 



3. Results 

3.1. Astrometric Confirmation 

Astrometric measurements are presented in Table [TJ First-epoch astrometry was measured 
from the companion centroid position relative to the starspot under the 0.3" mask in the "broadband 
ADI-only" reduction. The companion and starspot positions were measured in each frame, with 
the companion detected at S/N > 10 per frame. Because the companion falls on a slightly different 
part of the 0.3" mask in each image, angular asymmetries in the mask transmission can cause 
systematic astrometric trends of up to one pixel. Thus, we adopt conservative errors on our first- 
epoch astrometry based on the amplitudes of these trends. During acceptance testing, the NICI 
platescale was measured at 18 mas/pixel. 

For the second-epoch astrometry, we used the maskless ADI narrow-band dataset, measuring 
the centroid position of the companion relative to the unsaturated stellar peak in the final stacked 
image. The error was estimated from the astrometric rms scatter of the 3 simulated objects (§ 2.3) 
with Amag that best matched the observed object flux. Astrometric measurements from the 0.22" 
mask broadband and the 0.3" mask ASDI datasets yielded similar results. We adopt the direct 
ADI results to avoid uncertainties from the mask transmission calibration. 

In Fig. [21 we plot the measured astrometry as well as the expected motion if PZ Tel B were a 
background star at infinite distance, given the first epoch position and known proper and parallactic 
motion of the primary star. The motion of PZ Tel B over 13 months deviates from the background 
ephemeris at the 8.9a level. Astrometry from the 0.3" mask ASDI datasets for a fainter object 
~3.8" away from PZ Tel are also shown in Fig. [2j The relative motion for this other object changes 
as expected for a background object, confirming that published proper motion and parallax of the 
primary star combined with the astrometric calibration of NICI are accurate. 

We reduced archival VLT NACO data of this system ob tained on 22 July 2003 with the Ks- 
band filter and the 0.013"/pixel camera, originally reported in lMasciadri et al.l (|2005l ) . If the object 



detected by NICI was moving with respect to PZ Tel A as a background object, it should have had 
a separation of >0.5" in 2003 and would have been detected at >15cr in the VLT data. By scaling 
and adding a PSF star image (from an unsaturated November 2003 NACO dataset of Gl 86) to the 
NACO data, we estimate that PZ Tel B had a separation of <0.17" in 2003. 



3.2. Photometry, Mass Estimates, and Limits on Other Companions in the Field 



Aperture photometry was measured from both the 0.22" mask JHKg datasets and the direct 
ADI narrowband dataset using aperture sizes from 1-3 pixels in radius. Photometry results were 
consistent for all apertures; we adopt the 2-pixel aperture results (Tabled]). As described in Section 
2.3, simulated objects were inserted and retrieved to convert between measured flux and simulated 
Amag and also to determine photometric errors as a function of Amag. A JHKs color-color 
diagram for PZ Tel B compared to field objects is presented in Fig. [3l While NICI photometry for 
PZ Tel B is reported in the MKO filter system and photometry for PZ Tel A is reported in 2MASS 
magnitudes, the differences between the two systems are small compared to the measurement 
uncertainties. PZ Tel B's measured colors are similar to mid/late-M field dwarfs. 



W e estimate the mass and Tg// of the companion based on the DUSTY models of lChabrier et al, 
(J2OOOI ). PZ Tel B's J — H color is consistent with the DUSTY model colors expected from young 



substellar objects at 12 and 20 Myr, so use of t he DUSTY models is a ppropriate. Age estimates 

13±4 Myr for the /3 Pic 



for the /? Pic mo ving group include 12 4 Myr (IZuckerman et al 



2001 



member GJ 3305 (JFeigelson et al.ll2006l ). and 11.2 Myr (lOrtega et al.ll2009l ). To account for the age 
range cited in the literature, we estimated the mass using a uniform distribution of ages between 8 
and 20 Myr. We have chosen to work from the bolometric luminosity rather than absolute magni- 
tudes; bolometric lumin osities are less subscep tible to uncertainties in the model atmospheres than 
single band raagnitu des ( Chabrier et al.ll200d ). Comparing to colors of field M dwarfs compiled by 



Leggett et al.l (J20ld ). we estimate SpT=M7±2 for PZ Tel B. JHKs- Colors vary somewhat between 
field M dwarfs and younger objects; thus, we have placed conservative errors on our estimated spec- 
tral type. We estimate a bolometric magnitude for PZ Tel B of 10.96ib0.18 mag from our measured 
M-Kx and BC = 3.1ifc0.1 m ag (using the K-haxid bolometric correction vs. spectral type relation from 



Golimowski et al.ll2004al ). We derive a model mass of 36ib6 Mjup, model T^eff of 2702ib84 K and 



model log(g) of 4.20ib0.11 dex. Using the same input age grid and the single band magnitudes 
yielded similar mass estimates (see Figure [3l right p anel). Using the sam e bolometric magnitude 
methodology with the grainless NextCen models from lBaraffe et al.l (jl998l ). we derive similar values 
of 44ib9 Mjupand 2764ib67 K. Inter polating over the same age range, using our derived absolute 
//-band magnitude and the models of lBurrows et al.l (jl997l . I2OOII ) , we also find a similar mass range 
of 38±8 Mjup. 

Although a few fainter background objects were seen in the PZ Tel field at separations <8", 
no other common proper motion companion was detected by NICI. We achieved So" contrast limits 
of 12.9 mag at 0.5", 14.6 mag at 1", and 16.6 mag at 2.25" in our 1.6^m SDI -|- ADI dataset. These 
contrast l imits were verified by extensive simulated planet tests. Interpolating from the DUSTY 
models of lChabrier et al.l ()2000l ) and adopting an age of 12 Myr, our observations are deep enough 
to detect any planets with masses >6 Mj^p at separations >0.5". 



3.3. Constraints on Orbital Parameters 



We estimate the seminiajor axis of PZ Tel B's orbit from its observed separati on. Assuming 
a unif orm eccentricity distribution between < e < 1 and random viewing angles, iDupuv &: Liu 
([2010|) compute a median correction factor between projected separation and semimajor axis of 
1.10_g'3g (68.3% confidence limits). Using this, we derive a semimajor axis of 20_y AU for PZ Tel B 
based on its observed separation in May 2010. These correspond to an orbital period estimate of 



79_2i yr, for an assumed total system mass of 1.284; 



-0.050 
-0.200 



M,^. 



We used the 2009 and 2010 NICI astrometry for PZ Tel B to place constraints on its orbital 
eccentricity through Monte Carlo simulations that account for astrometric errors and all possible 
inclinations (0° < i < 90°). For each inclination, we used 10^ uniformly distributed values of the 
PA of the ascending node (0° < f] < 360°), allowing for different deprojections of the observed 
astrometry from the sky plane to the orbital plane. We computed the radial separation (r) and 
instantaneous velocity of PZ Tel B (v) for each possible orbital plane, using the parallax to convert 
to physical units. Assuming a system mass of 1.084-1.334 M©, we found the semimajor axis for 
each simulated orbit (a = r/(2 — ru^/Mtot))- By considering the amplitude of the transverse ve- 
locity of PZ Tel B in the orbital plane (vt), we found the eccentricity corresponding to each trial 

(e = A /l — r'^v^/a). In f«30% of cases, the trial orbits were unbound (e > 1, a < 0), and these were 
excluded from the analysis. Most of the resulting trial orbits are highly eccentric, and in Figured] 
we show the 2a minimum allowed eccentricity as a function of the assumed inclination. Our Monte 
Carlo simulations effectively rule out orbits with e < 0.6 for any assumed inclination, due to the 
fact that PZ Tel B displays a large radial velocity given the known mass of the system that cannot 
be totally accounted for by projection effects. We choose not to assume an input distribution for 
inclination that would allow us to derive confidence limits on eccentricity, as randomly oriented 
orbits are not likely to be appropriate because of observational selection effects. (Highly eccen- 
tric companions are preferentially dis covered face-on, wh ich is the opposite sense of the geometri- 
cal preference for edge-on orbits, e.g. iDupuy Sz Liull2010l ). 



The high eccentricity may be due to a dynamical interaction with another body. However, a 
radial velocity measurement bv ISoderblom et al.l (|l998l ) rules out the possibility that PZ Tel A is 



a dou ble-lined spectroscopic binary. Also, PZ Tel A shows a 70 //m emission excess (JRebull et al 



20081 ) , evidence for a significant quantity of circumstellar dust undisrupted by the orbi tal motion of 



the co mpanion. Assuming a simple blackbody and characte ristic temperature of 41 K.lRebull et al 



(|2008l ) set a minimum radius for the inner edge of 50 AU. lArtymowicz &: Lubowl (j 19941 ) find that 
binaries at a range of mass ratios and disk viscosities and with eccentricities > 0.6 will have -^^ > 
2.5 where re^^e is the inner edge of the circumbinary disk and a is the semimajor axis of the binary. 
If the inner edge of the circumbinary disk is at 50 AU, the semimajor axis of the orbit of PZ Tel 
AB is likely <20 AU. 



Conclusions 



With an estimated mass of 36ib6 Mj^p, PZ Tel B is a mong the lowest-mass com panions directly 
imaged around a young solar analog (e.g. compilation in lZuckerman &: Songil2009l ). There are only 
three other planetary or substellar co mpanions known to d ate in the ~10 Myr age range - the re- 
cently confirmed planet around /3 Pic ( Lagrange et al. 201Cll'l. the mid-to-late L dwarf 2MASS 1207- 



39B in the ;^8 Myr TW Hydra A ssociation (jChauvin et al.ll2005l ). and the ~20 Mjup late-M dwarf 



TWA 5B ([Lowrance et al.lll999l ). PZ Tel B will be a new benchmark companion for studies of 
ultracool, low-gravity photospheres. The projected separation of PZ Tel B is only 18 AU, making 
it one of very few substellar or planetary companions directly imaged at separations of <20 AU. 
Further astrometry and spectroscopy of this object will set additional limits on its orbital properties 
and provide improved estimates for effective temperature and surface gravity, better constraining 
the mass and formation history of this object. 
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Fig. 1. — Left: J, H, and Kg-hand images of the PZ Tel system obtained with NICI in direct 
imaging mode at the Gemini-South Telescope in May 2010. North is up, east is left. The primary 
resides at the center of the translucent 0.22" radius focal plane mask and is attenuated by a factor 
of 329 in J, 214 in H, and 131 in ii's-band. The confirmed companion is at 0.36" separation and 
PA=59.4° with flux ratios of AJ= 5.40±0.13, AH = 5.38±0.09, and AK^ = 5.04±0.10 mag. Right: 
Maskless ADI H2 2-1 (2.12 ^m) image from May 2010. Light from the primary has been removed 
by ADI processing. 



12 




2009 



2010 2011 

Epoch 

(a) PZ Tel B 



2012 




2010 201 

Epoch 



2012 



(b) Background Object 



Fig. 2. — On-sky motion measured for the PZ Tel B (a) and a faint background object at 3.8" 
in the same dataset (b). The expected motion and la uncertainty for a background object are 
overplotted, given the known proper and parahactic motion of the primary star, and the measured 
first-epoch position of the companion and background object. PZ Tel B cannot be a background 
object, as its motion diverges from the expected background track at the 8.9cr level. In contrast, 
the object at 3.8" separation moves as expected relative to the primary for a background object. 
Astrometry errors are measured in separation and PA relative to PZ Tel A and converted to RA 
and DEC. A 0.5° PA error at 3.8" separation translates to a considerably larger linear error than 
the same error at 0.33" separation, hence the larger error bars in RA for the background object. 
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Fig. 3. — Left: the JHKs colors of PZ Tel B compared to field objects ffrom lRavner et al 



200S 



figur e 37). The L dw a rf colors are dra^yn from a compilation by Sandy Leggett (ILeggett et al 
2nid : Ichhi et allbood : lOolimowski et al.ll2nn4bl : IXnapp etaLlbnoi lOeballe et allbood ). PZ Tel 
B's colors are plotted as a red circle and are consistent with those of a mid to late M dwarf. M 
dwarf points are shown as the a verage over seye r al ob jects at each subclass. Right: Mjj vs. age 
using DUSTY model values from IChabrier et al.l ([20001) • PZ Tel B is plotted as a red circle. 
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Fig. 4. — Monte Carlo constraints on inclination and eccentricity of PZ Tel B from our 2 epochs of 
NICI astrometry. For each inclination, 10^ values of the PA of ascending node (O) were randomly 
drawn to allow for different deprojections of the observed astrometry from the plane of the sky to 
the plane of the orbit. In r::30% of cases, the resulting trial orbits were unbound (e > 1, o < 0), 
and these were excluded from the analysis. Most of the resulting trial orbits are highly eccentric. 
These results constrain the orbit of PZ Tel B to eccentricities of >0.6. 
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Table 1. Properties of the PZ Tel AB System 



Primary 



Secondary 



Distance 

Age 

Proper Motion (/i„, ^s) 

Separation: 11 Apr 2009 UT 

Position Angle: 11 Apr 2009 UT 

Separation: 9 May 2010 UT 

Position Angle: 9 May 2010 UT 

A J (mag) 

AH (mag) 

AKs (mag) 

ACH4, 1% short (mag) 

AH2 1-0 (mag) 

J (mag) 

H (mag) 

Ks (mag) 

J — H (mag) 

H - Ks (mag) 

Mj (mag) 

Mh (mag) 

Mk, (mag) 

Spectral type 

Estimated Mass (from Lboi) 

Estimated Teff (from Lboi) 

Estimated log(g) (from hf,oi) 



51.5±2.6 pc"- 

12l^ Myr'' 

(17.6±1.1, -83.6±0.8) mas/yr'' 

0.330±0.010" (16.4±1.0 AU) 

59.0±1.0° 

0.360±0.003" (17.9±0.9 AU) 

59.4±0.5° 

5.40±0.14 
5.38±0.09 
5.04±0.15 
5.22±0.12 
5.12±0.13 
6.86±0.02^ 12.26±0.14 

6.49±0.05'^ 11.87±0.10 

6.38±0.02" 11.42±0.15 

0.37±0.05 0.39±0.17 

0.11±0.06 0.45±0.18 

3.30±0.12 8.70±0.18 

2.93±0.12 8.31±0.15 

2.82±0.12 7.86±0.19 

KO [M7±2]'^ 

Mo" 36±6 Mjup 
2702±84 K 
4.20±0.11 dex 



1.25 



+0.05 
'0.20 



van LeeuwenI (|2007f ) 



Zuckerman et al.l (|2001h 



'^from 2MASS 

"^Estimated from JHKs colors (§ 3.2) 



D'Antona fc Mazzitellil \l994 ) 



